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Abstract 

Along the line of thoughts of Berry and Robnik [1], the limiting gap distribu- 
tion function of classically integrable quantum systems is derived in the limit 
of infinitely many independent components. The limiting gap distribution 
function is characterized by a single monotonically increasing function fi(S) 
of the level spacing S, and the corresponding level spacing distribution is clas- 
sified into three cases: (i) Poissonian if /2(+oo) = 0, (ii) Poissonian for large 
S, but possibly not for small S if < p,(+oo) < 1, and (hi) sub-Poissonian 
if fi(+oo) = 1. This implies that even when the energy-level distributions 
of individual components are statistically independent, non- Poissonian level 

spacing distributions are possible. 
PACS number(s): 05.45.Mt, 03.65.Sq 
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I. INTRODUCTION 



An important property of quantum-classical correspondence appears in the statistical 
property of energy levels of bounded quantum systems in the semiclassical limit. Universal 
behaviors are found in the statistics of unfolded energy levels at a given interval [2-4], 
which are the sequence of numbers uniquely determined by the energy levels using the 
mean level density obtained from the Thomas- Fermi rule [5]. It is widely known that, for 
quantum systems whose classical counterparts are integrable (those systems will be referred 
to as classically integrable systems), the distribution of nearest-neighbor level spacing is 
characterized by the Poisson (exponential) distribution [3], while for quantum systems whose 
classical counterparts are strongly chaotic, the quantal level statistics are well characterized 
by the random matrix theory which gives level-spacing distribution obeying the Wigner 
distribution [5,6]. 

Level statistics for the integrable quantum systems has been theoretically studied by 
Berry- Tabor [3], Sinai [7], Molchanov [8], Minami [9], Bleher [10], Connors and Keating 
[11], and Marklof [12], and have been the subject of many numerical investigations. Still 
its mechanism is not well understood, the appearance of the Poisson distributions is now 
widely admitted as a universal phenomenon in generic integrable quantum systems. 

As suggested, e.g., by Hannay (see the discussion of [1]), one possible explanation would 
be as follows: For an integrable system of / degrees-of-freedom, almost every orbit is gener- 
ically confined in each inherent torus, and the whole region in the phase space is densely 
covered by invariant tori as suggested by the Liouville- Arnold theorem [13]. In other words, 
the phase space of the integrable system consists of infinitely many tori which have infinites- 
imal volumes in Liouville measure. Then, the energy level sequence of the whole system 
is a superposition of sub-sequences which are contributed from those regions. Therefore, if 
the mean level spacing of each independent subset is large, one would expect the Poisson 
distribution as a result of the law of small numbers [14]. This scenario suggested by Hannay 
is based on the theory proposed by Berry and Robnik [1]. 



2 



The Berry-Robnik theory relates the statistics of the energy level distribution to the 
phase-space geometry by assuming that the sequence of the energy spectrum is given by the 
superposition of statistically independent subspectra, which are contributed respectively 
from eigenfunctions localized onto the invariant regions in phase space. Formation of such 
independent subspectra is a consequence of the condensation of energy eigenfunctions on 
disjoint regions in the classical phase space and of the lack of mutual overlap between their 
eigenfunctions, and, thus, can only be expected in the semi-classical limit where the Planck 
constant tends to zero, H — > 0. This mechanism is sometimes referred to as the principle 
of uniform semi-classical condensation of eigenstates [15,16]. This principle states that 
the Wigner function of a semiclassical eigenstate is connected on a region in phase space 
explored by a typical trajectory of the classical dynamical system. In integrable systems, 
the phase space is folded into invariant tori, and the Wigner functions of the eigenstates 
tend to delta functions on these tori in the semiclassical limit [17]. On the other hand, in a 
strongly chaotic system, almost all trajectories cover the energy shell uniformly, and hence 
the Wigner functions of eigenstates are expected to become a delta function on the energy 
shell as suggested by the quantum ergodicity theorem [18,19]. Because of the suppression of 
the tunneling, each quantum eigenstate is folded into independent subsets in the semiclassical 
limit, and is expected to form independent spectral components. Indeed, formation of such 
independent components are checked numerically in a deep semi-classical regime [20]. 

In the Berry-Robnik approach [1] , the overall level spacing distribution is derived along a 
line of mathematical framework by Mehta [6] , as follows: Consider a system whose classical 
phase space is decomposed into iV-disjoint regions. The Liouville measures of these regions 
are denoted by pi(i = 1, 2, 3, • • • , N) which satisfy Yn=\Pi — 1- Let E{S) be the gap 
distribution function which stands for the probability that an interval (0, S) contains no 
level. E(S) is expressed by the level spacing distribution P(S) as follows; 



When the entire sequence of energy levels is a product of the statistically independent super- 





3 



position of N sub-sequences, E(S; N) is decomposed into those of sub-sequences, Ei(S; pi), 

N 

E(S;N) = Y[E i (S;p i ). (1.2) 

i=i 

In terms of the normalized level spacing distribution Pi(S] pi) of a sub-sequence, Ei(S; pi) is 
given by 

E i (S;p i )=p i da p t (x; pi)dx, (1.3) 

J S J (7 

and Pi(S] pi) is assumed to satisfy [1] 

r S- Pt (S; Pi )dS = -. (1.4) 
Jo Pi 

Equations (1.2) and (1.4) relate the level statistics in the semiclassical limit with the phase- 
space geometry. 

In most general cases, the level spacing distribution might be singular. In such a case, 
it is convenient to use its cumulative distribution function p,^ 

Pi(S) = / pi(x]pi)dx. (1.5) 
Jo 

In addition to equations (1.2) and (1.4), we assume two conditions for the statistical 
weights: 

• Assumption (i): The statistical weights of independent regions uniformly vanishes in 
the limit of infinitely many regions; 

maxpj^O as iV — > +oo. (1.6) 

i 

• Assumption (ii): The weighted mean of the cumulative distribution of energy spacing, 

Mp;A0 = f>/i,(p), (1.7) 



i=i 



converges in iV — > +oo to p(p) 



lim p(p;N)=p(p). (1.8) 

N— >+oo 



The limit is uniform on each closed interval: < p < S. 



In the Berry-Robnik theory, the statistical weights of individual components are related to 
the phase volumes of the corresponding invariant regions. This relation is satisfactory if 
the Thomas-Fermi rule for the individual phase space regions still holds, and thus supports 
eq.(1.4) [1]. Here we do not specify the validity of this problem and deal with the statistical 
weights as parameters. 

Under assumptions (i) and (ii), eqs.(1.2) and (1.4) lead to the overall level spacing dis- 
tribution whose gap distribution function is given by the following formula in the limit of 
N -> +00, 



where the convergence is in the sense of the weak limit. When the level spacing distributions 
of individual components are sparse enough, one may expect p — and the level spacing 
distribution of the whole energy sequence reduces to the Poisson distribution, 



In general, one may expect p 7^ which corresponds to a certain accumulation of the levels 
of individual components. 

In the following sections, the above statement is proved and the limiting level spacing 
distributions are classified into three classes. One of them is the Poisson distribution as 
discussed in the original work by Berry and Robnik [1] . The others are not Poissonian. We 
give examples of non-Poissonian limiting level spacing distributions in section III. In the 
concluding section, we discuss some relations between our results and other related works. 




(1.9) 



Ep^S) = exp (-S). 



(1.10) 



II. LIMITING LEVEL SPACING DISTRIBUTION 



A. Derivation of the limiting Gap distribution 



In this section, starting from eqs.(1.2) and (1.4), and assumptions (i) and (ii) introduced 
in the previous section, we show that, in the limit of infinitely many components N — > +00, 
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the gap distribution E(S; N) converges to the distribution function (1.9) with Jx. The 
convergence is shown as follows. 

Following the procedure by Mehta(see appendix A. 2 of Ref. [6]), we rewrite E(S; N) in 
terms of the cumulative level spacing distribution Pi(S) of independent components: 



N 



E(S;N) = H 



i=i 



da {I - fii(a)) 



(2.1) 



Equation (1.4) and integration by parts lead to 

r+oo r+oo | 

/ da{\ - mio)) = o{\ - mio))^ 00 + ap t (a)da = -, 
Jo Jo Pi 

where lim^+oo <r(l — Pi(a)) = follows from the existance of average, and hence to 

r+oo rS 

Pi da(l - pi(<j)) = 1 - pi / da(l - 
Js Jo 



(2.2) 



Since the convergence of J2f=i PiPi(cr) — > p(a) for N — > +oo is uniform on each interval 
a G [0, S] by Assumption (ii), E(S; N) has the following limit, 

N 

log£0S;A0=£log 



= -E 



1- Pi [ da(l - Hi(a)) 
Jo 

rS 

Pi / da(l - p l {a)) + 0(p 2 ) 
Jo 

N 

- / Ar[i-p(a ; A0] + £o(p 2 ) 
Jo T 

— > — / dcr [1 — /i(cr)l as iV — > +oo, 

JO 



(2.3) 
(2.4) 



where we have used \pi{u)\ < 1, log(l + e) = e + 0(e 2 ) in e <C 1, and the following property 
obtained from Assumption (i), 

N N 

| 0(p 2 ) I < C • max pj • p^ = C ■ max p^ — > as iV — > +oo, (2.5) 
i=i «=i 

with C a positive constant. Therefore, we have the desired result: 

rS 



lim E(S;N) = EJS) = exp 

N— »+oo 



[ (l-p(a))da 
Jo 



(2.6) 
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B. Weak convergence limit of the level spacing distribution 



In this section, we show that in N — > +00 limit, the level spacing distribution P(S; N) 
converges weakly to Pp(S), 



d 2 



(2.7) 



dS 2 MV 

According to the Helly's theorem [21,14], the weak convergence of the level spacing 
distribution is defined by 

rS r S 



^lim J P(x;N)dx = J Pp,(x)dx. 



(2.8) 



Since each side of the above equation is rewritten as, 

1 s 



( P(x;N)dx = 
Jo 

the limit (2.8) is equivalent to 



J 



s 



P n (x)dx 



d . « 



The above equation is proved as follows: By using equation (2.1), we rewrite -^E(S; N) in 
terms of the cumulative level spacing distribution function (j,i(S) of spectral components, 



d 

Is 



N 



1 - fjti(S) 



E(S;N) = -E(S;N)J2- , 

i=i Pi Is da{l - (m((t)) 



In the limit N — > +00, one has E(S; N) — > E^S) as shown by equation (2.6) and, 



N 



1 Pi Is °° da(l - Hi(a)) i=1 1 - Pi J da(l - 



N 



N 



i-^pMS) + ^o(p1) 



i=l 



1=1 



— ► 1 — fi(S) as N — > +00, 

where we have used equation (2.2), 1/(1 — e) = 1 + 0(e) in e < 1, and the limit (2.5). 
Therefore, we have the desired result: 



Jsz. £wo = -(i- m) E f ( S ) = ±e>(8). 



(2.9) 



We remark that, when the limiting distribution function p(S) is different iable, the asymp- 
totic level spacing distribution is described as follows: 

rS 



exp 









[1 — /x(cr)) da 



(2.10) 



C. Properties of the limiting level spacing distribution 

Since fJ,i(S) is monotonically increasing and < fJ>i(S) < 1, fi(S) has the same properties. 
Then, 1 — p(S) > for any S > and one has 
1 r s 

— da(l - p,(a)) — >l-fi(+oo) as S -> +oo. (2.11) 
S Jo 

According to the above limit, the level spacing distribution is classified into the following 
three cases in the sense of weak limit: 

• Case 1, p,(+oo) = 0: The limiting level spacing distribution is the Poisson distribution. 
Note that this condition is equivalent to p,(S) = for V S because Jx(S) is monotonically 
increasing. 

• Case 2, < p,(+oo) < 1: For large S values, the limiting level spacing distribution 
is well approximated by the Poisson distribution, while, for small S values, it may 
deviate from the Poisson distribution. 

• Case 3, p,(+oo) = 1: The limiting level spacing distribution deviates from the Pois- 
son distribution for V S, and decays as S — > +oo more slowly than does the Poisson 
distribution. This case will be referred to as a sub-Poisson distribution. 

One has Case 1 if the individual cumulative distribution function ^(S) are bounded by a 
finite positive function g(S), 

ia(S) < pUS), 

for y i and r\ > 0. Indeed, one has, 

N N 

f,(S;N) = < g(S)Y,PiPi < — = ft(S), as N oo. 

i=i % 
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More specifically, for example, one has Case f if the individual level spacing distributions 
are derived from scaled distribution functions fa as 1 , 

Pi(S;pi) = pifi(piS), 

where /j satisfy 

r+oc r+oo 

/ fi(x)dx = 1, / xfi(x)dx = 1, 
Jo Jo 

and are uniformly bounded by a positive constant D: \fi(S)\ < D (1 < i < N and S > 0). 
Indeed, one then has 



N fS 

HS;N)\<- 



Pi / Pi(%,Pi)dx 

i=l J ° 

N g 

<Y.Pi \fi(pix)\dx 
i=i Jo 



(2.12) 



N N 

< DSJ2p1 < DSmaxpiJ^Pi — > = ft(S). (2.13) 

i=i * i=i 

This includes the case studied by Berry and Robnik [1], where the gap distribution 
is a product of superposition of a single regular component characterized by the Poisson 
distribution and iV equivalent chaotic components characterized by the Wigner distribution, 
and the latter is expressed by the product of the scaled distributions as: 

N 

E BR (S; N) = exp (-p S) JJ Sf IGNER (5; Pi ), (2.14) 

i=i 

where the statistical weights are pi = ^-jp- and the individual level spacing distributions 
fi corresponding to the gap distributions £' J WIGNER (S'; Pi ) = erfc {^- Pi S^ are given by the 
dimensionless Wigner distribution: 



1 Eq.(1.4) is described by rewriting x = piS in the following way, 1 = / SpiPi(S; pi)dS = 
j xfi(x)dx, so that fi(x) = PiPi(x; pi) is a dimensionless function. For instance, fi(x) = exp(— x) 
in the case of Poisson distribution pi(S; pi) = pi exp (—piS), and f(x) = ^xexp (—jx 2 ) in the case 
of the Wigner distribution pi(S; pi) = ^pjSexp [-^p 2 S 2 ] . 
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fi( x ) = -y exp 



71 2 



Indeed, one has the Poisson distribution in iV — > +00 limit: 

E BR {S; A) = exp 



- Po S + (N - 1) logerfc [ ^\^S 



2 A- 1 



(2.15) 



(2.16) 



III. EXAMPLE 

As an example of the deviation from the Poisson distribution, we study rectangular 
billiard system whose energy levels are described by using positive integer numbers, m and 
% as follows, 



2 *2 



(3.1) 



where a is the aspect ratio of two sides of a billiard wall. For a given energy interval [e, e+ Ae] , 
each energy level is classified into components according to the eigenvalues, % — 1, 2, 3, • • • , N; 



N = 



'1 + 7 



a 



(3.2) 



where 7 = Ae/e, and [x] stands for the maximum integer which does not exceed x. The 
relative weight of each component is given by 

fs-(*) 2 )+°0 



Pi 



4(1+7) 



1 - 



JV 



if i < 



N 
/T+7 ; 



4(1+7) 
Nttj 



(3.3) 



if 



JV 

/ T+7 



< i < N. 



As easily seen, satisfies the assumption (i); 







as A — > +00. 



(3.4) 



Note that the limit of infinitely many components, A — > +00, corresponds to the high 
energy limit, e — > +oo(see Eq.(3.2)), which is equivalent to the semiclassical limit in physical 
systems. In this limit, the statistical weight of each sub-spectrium becomes sparse, since 
each element of /i(S] A), pi/i^S), tends to zero: piPi(S) < max.,- pj — > 0. 
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Figures 1(a) and 1(b) show numerical results of the level-spacing distribution P(S) for 
two values of a, and figures 2(a) and 2(b) show the gap distribution function corresponding 
to figures 1(a) and 1(b), respectively. In case that a is far from rational, P(S) and E(S) 
are well approximated by the Poisson distribution while in case that a is close to a rational 
expressed as a = p/q, where p and q are coprime positive integers, they deviate from the 
Poisson distribution [22]. 

In order to compare the non-Poisson distribution and the classification given in the 
previous section, we consider, 

^"- ''Kr ' (3 - 5) 

When iV — » +00, p,(S;N) approaches p,(S), and this function distinguishes the three 
cases as follows: In Case 1 i.e., where the level spacing obeys the Poisson distribution, 
limjv^ +00 p,(S; N) = 0. In Case 2, \im N ^ +OQ £l(S; N) — > c as S — > +00 (0 < c < 1), and in 
Case 3, where the sub-Poisson distribution is expected, \im.N_> +00 £l(S; N) — > 1 as S — > +00. 

Figure 3 shows p,(S; N) for different values of N. The dotted line p, — exhibits the 
Poisson distribution. From this, one can think that jj(S;N) for iV = 61905,5* < 10 well 
approximates \im.N^ +OQ £l(S; N). 

Figure 4 shows jj,(S;N) for the two values of a corresponding to figures 2(a)-2(b), re- 
spectively. In case that the numerical data is well characterized by the Poisson distribu- 
tion (figures 1(a) and 2(a)), the corresponding function p,(S;N) agrees with 0, while in 
case that deviates from the Poisson distribution (figures 1(b) and 2(b)), jj(S;N) ^ and 
jj,(S; N) — > c(0 < c < 1) for S — > +00. Therefore, this result corresponds to the Case 2. 

In this model, we have not yet observed the clear evidence of Case 3. Such a case is 
expected when there is stronger accumulation of the energy levels of individual components. 

IV. EXTENDED FORMALISM OF THE BERRY-ROBNIK DISTRIBUTION 



In this section, we propose one possible extension of the Berry- Robnik distribution (2.14) 
for the level statistics of the nearly-integrable system with two degree-of-freedom. Since the 
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classical phase space of the nearly-integrable system consists of regular and chaotic regions 
and the Liouville measures of the chaotic regions are larger than zero, pi > 0, this system 
does not support the assumption (i). However, the regular regions consist of infinitely 
many subsets, and our approach shown in section II is partially applicable to the spectral 
components corresponding to the regular regions. 

Following to the assumption proposed by Berry and Robnik [1], the gap distribution 
functions in the nearly integrable system, which are contributed from the individual chaotic 
regions, are characterized by the Random Matrix Theory(RMT) . Then one has, 

JVi N 2 

\ogE(S-N 1] N 2 )=\ogl[Ef M ^S;p l )l[E J (S;p J ) (4.1) 

i=l j=l 
Ni N 2 

= J2\ogE?^(S;p t )+j:\ogE J (S;p J ), (4.2) 

i=l 3=1 

where p'j = Pj+N^ and Ej(S; p'j) denote the gap distribution functions corresponding to the 
subsets in the regular regions. As shown by eq.(2.6), Ej has the following limit, 

N 2 , rS 

Y 

3 = 

where p = J2f=i Pj and p(cr) = \im N2 ^ +O0 J2j ! =iPjPj( <J )- Accordingly, in the partial limit of 
Ni <C +oo, N 2 — > +oo, the original proposal for the gap distribution by Berry and Robnik 
is replaced by 

rS 



N 2 S 

log Ej(S; pj) — > p / da [I- p(a)\ as N 2 -> +oo, 



Jim E(S; iV i; N 2 ) = E^S; N,) = exp 



-po ( (1 - p(a))da 
Jo 



Ni 

UEt MT (S;Pi), (4-3) 

J 1=1 



where p denotes the total amount of the Liouville measures of the regular region in mixed 
phase space. The above distribution formula is classified into the following three cases; 
Case 1', p(+oo) = 0: Berry-Robnik distribution, Case 2', < p(+oo) < 1: Berry-Robnik 
distribution for large S, but possibly not for small S, and Case 3', p,(+oo) — 1: A distribution 
function obtained by the superposition of spectral components obeying the sub-Poisson 
statistics and the Random matrix theory. From this classification, one can see that the new 
formula (4.3) admits deviations from the Berry-Robnik distribution when p(+oo) ^ 0. 
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We remark that P(S; N 1: N 2 ) = -£zE(S; N 1: N 2 ) in the limit N 2 -> +00 converges weakly 
to the limiting level spacing distribution: Pp,(S; N±) = -^E^S; Ni), and when the limiting 
function fi(S) is different iable, the asymptotic level spacing distribution admits the following 
density: 



exp (-ft j\l -fl(a))d^j n^f MT (5; Pl ) 



(4.4) 



dS 2 

The validity of the Berry-Robnik distribution has been checked for generic nearly- 
integrable systems by many numerical investigations [23-29]. Among them, Prosen and 
Robnik found numerically for several systems that there is a high energy region in which the 
Berry-Robnik distribution formula (2.14) well approximates the level spacing distribution 
[23]. This energy region is sometimes referred to as the Berry-Robnik regime [20]. While 
they also found that the level spacing distribution in the low energy region deviates from 
the Berry-Robnik formula, and approximates the Brody distribution. This behavior was 
studied in terms of a fractional power dependence of the spacing distribution near the origin 
at S = 0, which could be attributed to the localization properties of eigenstates on chaotic 
components [23,24]. From the above classification, Case V:p,(+oo) = sould be satisfied in 
the Berry-Robnik regime. While Case 2' and Case 3' might propose another possibilities. 
When the spectral components corresponding to regular regions show strong accumulation, 
the gap distribution obeys the distribution formula (4.3) with < /2(+oo) < 1, and the level 
spacing distribution shows deviations from the Berry-Robnik distribution. Therefore, this 
result might propose another possibility of the Berry-Robnik approach. 



V. CONCLUSION AND DISCUSSION 

In this paper, we investigated the gap distribution function of systems with infinitely 
many independent components and discussed the level-spacing statistics of classically inte- 
grable quantum systems. In the semiclassical limit, reflecting infinitely fine classical phase 
space structures, individual energy eigenfunctions are expected to be well localized in the 
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phase space and contribute independently to the level statistics. Keeping this expectation in 
mind, we considered a situation in which the system consists of infinitely many components 
and each of them gives an infinitesimal contribution. And by applying the arguments of 
Mehta, and Berry and Robnik, the limiting level spacing distribution was obtained whose 
gap distribution function is described by a single monotonically increasing function p,(S) of 
the level spacing S: 



The weak convergence limit of the level spacing distribution is classified into three cases; 
Case 1: Poissonian if p,(+oo) = 0, Case 2: Poissonian for large S, but possibly not for small 
S if < p,(+oo) < 1, and Case 3: sub-Poissonian if /2(+oo) = I. Thus, even when the 
energy levels of individual components are statistically independent, non-Poissonian level 
spacing distributions are possible. 

In most general cases, the integral in equation (5.1) converges in S <C +oo and then 
linis^+oo Efi(S) 7^ 0, the limiting gap distribution Ep(S) does not work accurately. In such 
case, however, its differentiation (2.9) still work accurately in S — > +oo limit [30], and thus 
the above classification (Case 1-3) holds in general. 

Note that the singular level spacing distribution can be taken into account in terms of 
non-smooth cumulative distributions. Such a singularity is expected when there is strong 
accumulation of the energy levels of individual components. For a certain class of systems, 
such accumulation is observable. One example is shown in section III where the results 
show clear evidence of Case 2. Another example is the two-dimensional harmonic oscillator 
whose level spacing distribution is non-smooth for arbitrary system parameter [3,10]. The 
final example is studied by Shnirelman [31], Chirikov and Shepelyansky [32], and Frahm and 
Shepelyansky [33] for a certain type of system which contains a quasi- degeneracy result from 
inherent symmetry(time reversibility). As is well known, the existence of quasi-degeneracy 
leads to the sharp Shnirelman peak at small spacings. 

Finally, in section IV, we proposed one possible extension of the Berry-Robnik distribu- 




(5.1) 
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tion for classically nearly-integrable quantum systems. This extension admitted deviations 
from the Berry-Robnik distribution when there is strong accumulation of the energy levels 
of spectral components. Such possibilities will be studied elsewhere. 
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FIGURES 
(a) (b) 




S S 

FIG.l Numerical results of the level spacing distribution P(S) for (a) a — 1 + | x 10~ 4 , 
(b) a = 1 + f x KT 9 . We used energy levels e m ,i with ^e™,, G [300 x 10 7 , 301 x 10 7 ]. Total 
numbers of levels are (a) 10000016, (b) 10000046. The dotted curve in each figure shows 
the Poisson distribution: P(S) = e~ s . 
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FIG.2 The gap distribution function E(S) for (a) a = 1 + f x 10~ 4 , (b) a = 1 + f x 10~ 9 . 
The dotted curve in each figure exhibits the Poisson distribution: E(S) = e~ s . 
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FIG.3 jl(S; N) for N = 2037 and for N = 61905. In each case, we fixed a = 1 + f x 10~ 4 , 
and used 30048 energy levels with ^e mtl G [323 x 10 4 ,326 x 10 4 ] and 10000016 energy 
levels with 4^e m ,« £ [300 x 10 7 , 301 x 10 7 ], respectively. The dashed line (p, = 0) exhibits 
the Poisson distribution. 
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FIG.4 The distribution function p(S; N) for (a) a = 1 + f x 10~ 4 (N = 61905), and (b) 
a = 1 + f x 10~ 9 (iV = 61906). The dashed line (/2 = 0) exhibits the Poisson distribution. 
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